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; ABSTRACT 

Qh' We discuss a semi-analytical solution of the transport equation for electrons at 

O ' a non-relativistic shock in the presence of synchrotron energy losses. We cal- 

C/3 , culate the spectrum of accelerated (test) particles at any point upstream and 

downstream of the shock for an arbitrary diffusion coefficient and we specialize 
the results to three cases: 1) diffusion constant in momentum {D{p) — Dq), 
' 2) Bohm diffusion {D{p) cx p), and 3) Kolmogorov diffusion {D{p) cx p^^^). 

I Of special importance is the determination of the shape of the cutoff in the 



electron spectrum which depends on the diffusion properties felt by particles 
in the shock region. The formalism can be generalized to the case of a shock 
with an upstream precursor induced by the dynamical reaction of accelerated 
particles. 

Key words: acceleration of particles - shock waves 



1 INTRODUCTION 

Particle acceleration of charged particles in the presence of losses has been subject of much 
investigation in the past, but in recent times its importance has been revived by the detection 
of non-thermal radiation from supernova remnants (SNRs) that, directly or indirectly, can be 
attributed to accelerated electrons in strong magnetic fields or intense radiation fields. Two 
instances can help stressing this point: observations carried out with the Chandra, XMM 
and Suzaku telescopes have shown that the non-thermal X-rays from SNRs are the result of 
synchrotron emission of highly relativistic electrons in magnetic fields which most likely are 
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amplified with respe c t to the interstellar fie 



Lazendic et al. 



2003 



2004 



Vink fc Laming 



ds by a factor ^ 1 00 (IBamba et al. 



2003 



2003 



Ballet 



2005 : 



20061 ) ■ In some cases the spectra 



of this radiation show a pronounced cutoff corresponding to the maximum energy of the 
radiating electrons. The spectrum of the X-rays and especially the shape of the cutoff carries 
precious information on the acceleration process and on the possible shock modification 
induced by the dynamical reaction of accelerated particles. Moreover, the morphology of the 
X-ray emission may lead to a direct estimate of the magnetic field strength in the shock 
region, which in turn affects particle diffusion. 

The magnetic field strength and topology determines the p ossibility for hadronic cosmic 



rays to reach the knee energy (IBlasi. Amato fc CaprioU 



20071 ). In turn the ma gnetic field 



1978 . 



can b e generated by high energy cosmic rays if they excite streaming instabilities (iBelll . 
2OO4I ). Electrons are affected by the magnetic field through diffusion and energy losses: both 



these effects may produce visible features in the cutoff region of the electron spectrum and 
eventually manifest themselves in the spectrum of the emitted radiation. 

Here we discuss the process of test-particle acceleration at a non-relativistic shock in 
the presence of synchrotron losses for three choices of the diffusion coefficient: 1) diffusion 
constant in momentum {D{p) = Dq), 2) Bohm diffusion {D{p) oc p), and 3) Kolmogorov 
diffusion {D{p) oc p^^^). The transport equation at the shock surface is transformed in a 
simple equation to be solved iteratively and returns the spectrum of accelerated particles 
at the shock, while the spectrum at any location upstream and downstream is determined 
analytically through the proper Green functions. 

Clearly this is not the first calculation of the process of electron acceleration at shocks, 
but previous methods that one can find in the literature are either numerical (solution of 
the transport equation through finite differences schemes) or analytical but applied to spe- 
cific choices of the diffusion coeffici ent (D(p) constant in moni e ntum) or to limited regions 
of momentum space. For instance IWebb. Drury fc BiermannI (119841 ) used a Green func- 



tion method to solve ana^ 
was then investigated by 



ytically the ca s e of c onstant diffusion coefficient. The same case 



Webb fc Fritzl (Il990l ) in the context of th e so-ca lled generations 



method, an elegant generalizatioii of the stochastic method of iBelll ( 119781 ) for accelerated 



protons. 



Zirakashvili &: AharonianI (120071 ) found an analytical solution for the case of Bohm 



diffusion but limited to particle momentum larger than the maximum momentum. The 
authors proposed a fitting formula, derived however from a numerical (finite differences) 
calculation, in order to interpolate the low energy and high energy regimes. The work by 
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Zirakashvili fc AharonianI (120071 ) is particularly important since it emphasizes the fact that 



in the case of Bohm diffusion, the spectrum of electrons has a super-exponential cutoff. This 
has obvious consequences in the comparison with observations of X-rays due to synchrotron 
emission of high energy electrons in the cutoff region. A numerical investigation of the prob- 
lem, also including inverse Compton losses in the Klein-Nishina regime, was presented by 



Vannoni. Gabici fc AharonianI (120091 ). In addition to these approaches that aimed at solving 



the exact transport equation, ther e are also attempt s to face the problem of ele ctron accel 



eration within leaky-box models (jPrury et al 



1999 



Protheroe fc Stanev 



19991 ). A generic 



prediction of these models is the appearance of pronounced pile-ups in the spectra of accel- 
erated particles when the compression factor is larger than 4. Such pronounced and sharp 
features are typically absent in more detailed calculations. 

The method discussed here provides both the spectrum of accelerated particles and their 
spatial distribution upstream and downstream of the shock, which is important to determine 
the morphology of the synchrotron emission. The calculation can in principle be carried out 
for any diffusion coefficient as function of momentum. 

In this paper we restrict ourselves to the test particle regime of electron acceleration at 
a non-relativistic shock. However, while doing so we will obtain some results that will later 
be used to generalize the calculations to the case of cosmic ray modified shocks. 

The paper is structured as follows: in [|2] we illustrate the general theoretical bases of 
the method and the procedure to solve the relevant equations. In ^ we discuss the imple- 
mentation of the iterative procedure, and illustrate the results of the calculations in terms 
of spectrum of the accelerated particles and their spatial distribution around the shock. We 
conclude in © 



2 THE TRANSPORT EQUATION AND GREEN FUNCTIONS 



The acceleration of electrons at a shock front is described by the following transport equation: 



u 



df 



dx 
where p 



d_ 
dx 



Dip) 



df{x,p) 
dx 



1 dzL d fix n) 1 d 

+ ol~P ^ 2^ [PP^fi^^P)] + Q{x,p), 

6 ax op p"^ op 



(1) 



—Ap^ {A = ^arU B / {rrtecY , with Ub = -B^/Svr, B the magnetic field, ctt the 



Thomson cross section, mg the electron mass) is the rate of momentum loss due to syn- 
chrotron emission, f{x,p) is the distribution function of accelerated particles and Q is the 
injection term, assumed to have the form: 
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Q{x,p) = -^^S{p - p')S{x - x'). 



(2) 



In Eq. [T]we assume that energy losses allow the system to reach stationarity {df/dt = 0). 
In the test-particle regime, there is no precursor, therefore as usual du/dx = {u2 — Ui)5{x), 
and the x axis is oriented from upstream (— oo ^ a; ^ 0) to downstream (0 ^ x ^ +oo). 

For our purposes it is useful to introduce the function N{x,p) = Anp^ f{x,p), so that the 
transport equation reads: 



d_ 

dx 



uN - D{p) 



dN 
dx 



d_ 
dp 



[Ap'N] +Qo6{p-p')6{x-x'). 



(3) 



The solution can be found by requiring the continuity of the distribution function and 
by satisfying the boundary condition at the shock, which is easily found by integrating the 
full transport equation in a narrow region around the shock. We will return to this point 
later. 

We now focus our attention on solving Eq. [3]in the d ownstream and upstream fluids sep- 
arately. For this purpose, following Webb fc Frit2 f 199ol ). we multiply Eq. [3] by an arbitrary 
function Q{x,p; x',p') and integrate both terms for xi ^ x ^ X2 and pi ^ p ^ P2, where the 
extremes of integration will be defined later in a suitable way. This results in the following 



expression: 

rP2 



dp 



pi 



ON 

QiuN-D-) 



X2 



+ 



P2 



dp 



pi 



Xl 



X2 Q 

dxD — 
ox 



dx 



r dx [Ap'Ng]^] 

Xl 



P2 



dp 



pi 



X2 



Xl 



dxN i 

ox 



D 



dx'^ 



2dQ 



op 



0, 



(4) 



where the injection term is not present as long as we concentrate on the regions upstream 
and downstream of the shock and we assume that the injection occurs exactly at the shock 
surface. 

Here we also made use of the identity 



dNdg _ d 

dx dx dx 



dx 



dx"^ 



(5) 



We recall now that the function g{x,p] x' ,p') was chosen to be an arbitrary function of 
its arguments. We have therefore the freedom to choose it in the way that is most useful for 
our purposes. We choose the function g as the solution of the equation 



+ - Ap'^ = -6{x - x')6{p - p') 



2dG 



dx 



dx^ 



(6) 
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In the downstream fluid we choose Xi = and X2 = +00, while in the upstream fluid, 
Xi = — 00 and X2 = 0. In both cases pi = and p2 = 00. Eq. [6] is sometimes called the 

adjoint equation. Using these conditions in Eq. |4]we obtain the two fundamental expressions: 

dQi{x',p';x,p) 



N,{x,p) = - / dp'D^{p')No{p') 
Jp 

in the upstream region, and 



dx' 



N2{x,p) 



dp'D2{p')Noip') 



dQ2{x',p';x,p) 
dx' 



+ 



dp'Qi (x, p; X, p) ij>') (7) 



dp'g2{x',p;x,p)\^^^o(f)2ip'), (8) 



in the downstream region. In order to keep the notation as simple as possible we inverted 
all primed and unprimed variables, so that the physical quantities are all expressed as func- 
tions of unprimed variables. The indexes "1" and "2" refer respectively to quantities in 
the upstream and downstream fluid. Nq{p) is the spectrum of accelerated particles at the 
shock location. Clearly the distribution function of the accelerated particles is continuous 
across the shock, namely lim^^Q- Ni{x,p) = limx_,Q+ N2{x,p) = No{p). The lower limit of 
integration in Eqs. [7|and[H]is p instead of zero because of the definition of the function Q, 
solution of Eq. [6] (see Appendix). From the physics point of view, it is obvious that it must 
be so, because the contribution to the spectrum at a given momentum p can only come from 
particles with larger momentum that lose energy through synchrotron losses. The functions 
are defined as: 



01,2 (p) = ui^2No{p) - £'1,2 (p) 



dN{x,p) 
dx 



1,2; 



(9) 



which only depend on quantities evaluated at the shock. 

It should be noted that Eqs. [7] and [8] define the spectrum at any location upstream and 
downstream in terms of No{p) and once the functions Qi and Q2 are known. The solution 
at the shock, Nq{p), must be derived by using the boundary condition at the shock location. 
On the other hand, the Green function of the adjoint equation, Q, is easily calculated (see 
Appendix) to be: 



Q{x',z';x,z) 



TT 



2nAy t{z', z) 



exp 



{x' -x + {u/A){z- z')y 
4t{z', z) 



Here we introduced z = 1/p and the function r, which is defined as follows: 



T[Z , Z] 



r 



(10) 



(11) 



The main difference between the approach described here an d the approach of 



Webb. Drurv fc B' 



(jl984j ) is in the fact that our approach (previously mastered by lWebb Sz Fritzl (ll99ol ) for dif- 
fusion coefficient independent of momentum) allows one to write the distribution function 
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N{x, t) as an integral on momentum of the distribution function at the shock, No{p), and the 
particle flux weighed with the appropriate Green functions, which are easy to calculate. 
This proves convenient from the point of view of adopting iterative techniques to solve the 
equations. 

In Eq. [To] all dependence on the diffusion coefficient is contained in the function r, which 
can be calculated for any choice of D{jp). The Green function calculated here is the one 
which vanishes at upstream and downstream infinity, but not at the shock. One could have 
chosen the Green function that satisfies the boundary condition of vanishing at the shock 
surface and the equations for N{x,p) would have been simpler. However the calculations of 
such Green function is simple and co mputationally conve nient only for the case of constant 



diffusion coefficient (see for instance (IWebb fc Fritz 



1990)). 



The boundary condition at the shock is easily obtained by integration of the full transport 
equation around the shock surface, and reads: 



A(p) 



dN 
dx 



Di{p) 



dN 
dx 



1 dN^ 

+ -{U2 - Ui)p-— 

6 dp 



'iu2-ui)No = 0. 



(12) 



We introduce now the two functions 



Up) 



Dijp) dN^ 
uiNo{p) dx 



\x->Q- 



and 



D2{P) 



ID ( 1 

dpD^ip)No{p)n,--^lim 

D^ip) dN 



dpNo{p)Ci [1-Xi 



(13) 



lim 



U2No{p) x^O 

where 

7^1,2(X, Z, Z') 
Cl^2{x,Z,z') 



Up) 

dpD2{p)Noip)n2 + 



U2Nq{p) dx 

D2{P) 



lim 



d_ 

dx 



dg 

dx 



No{p) -^ojp 
QUx = 0) 



2r 



dGUx = 0) _ gi,2{x = 0) 



dpNo{p)C2[l-Up)] 

{-X+{U,,2/A^,2){Z-Z'))' 



2t 



dx 



2r 



, '"1,2 / 

-x + -—{z 

^1,2 



(14) 

(15) 

(16) 

0) one 



It is important to notice that in the limit in which losses are negligible (r 
has Xi ^ 1 and I2 —>■ 0. This behavior is exactly what one would expect based on the 
standard test-particle solution of the transport equation in the absence of losses: a physically 
meaningful, stationary solution can in fact be obtained only if the distribution function is 
homogeneous downstream, namely = 0. At the same time the boundary condition 



upstream of the shock reads Di 



dN{x=0-) 
dx 



-uiNq. This is exactly what the equations above 
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are confirming. On the other hand, when synchrotron losses become important and the limit 
r ^ does not apply, Eqs. [13] and [TH allow us to write the natural generalization of the 
jump conditions to the case with synchrotron losses. 

In terms of Xi 2, implicit functions of Nq{p), one can rewrite the boundary condition 
equation in a way which is reminiscent of the standard boundary condition for the loss-less 
stationary test-particle case: 

p y \{u2-ui) 

In the loss-less case, X2 — and Xi ^ 1, so that Eq. [17] reduces to the equation for the 
slope of a power law, which would then be s = (2^2 + Ui)/{u2 — Ui) = —(2 + r)/(r — 1) 
where r = Ui/u2 is the compression factor at the shock. The apexes k and k — 1 which have 
been introduced in Eq. [T7]are motivated by the fact that we intend to solve this equation by 
iterations: at each iteration the functions Xj are evaluated at the previous step. The method 
converges in a few iterations to the spectrum of particles at the shock location, Nq{p). At 
that point, the distribution function at an arbitrary point upstream and downstream is fully 
determined through Eqs. [7] and [3 The method is computationally much faster than any 
finite differences scheme of numerical integration of the transport equation. 




3 MAXIMUM MOMENTUM 



In the presence of losses, the maximum momentum can be estimated by requiring equality 
of the acceleration time and the loss time, the latter weighed by the propagation times in the 
upstream and downstream regions. The condition of balance between losses and acceleration 
reads: 

'dE\ 



'dE\ 

/ loss,! 



-TT - r« 2 = 



dt) 



(18) 



loss,2 



where tr are the residence times in the upstream and downstream regions. If we use, following 



torury 



19831), r_ = E/ (f)^ 



UI—U2 



Di{p) I D2{p) 

Ul U2 



condition 



Taccip) 



TR,1 + TR,2 
TR,! _|_ ''"i{,2 ' 



as the acceleration time, one has the 



(19) 



where tioss = E/{dE/dt)ioss is the loss time in the upstream (1) and downstream (2) re- 
gion. This condition defines the maximum momentum Pmax ioT any choice of the diffusion 
coefficient. 
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In the following we provide some useful expressions for the three choices of the diffusion 
coefficient used in the present calculation. Some more details are provided for the case of 
Bohm diffusion since it is often adopted for phenomenological applications, especially when 
the magnetic field is strongly turbulent. 

3.1 The case of constant diffusion coefficient 

For a constant diffusion coefficient D{p) = Dq, we have 

t{z,z')^^{z-z'). (20) 

In principle, the diffusion coefficient and the rate of losses can be different upstream and 
downstream. This generalization is trivial to implement if needed and reflects in using the 
appropriate constant values for Do and A. 

The maximum momentum is defined by the expression: 

p - ' "'^^-^^ (21) 
6DQAr r + 1 

where r — Ui/u2 is the compression factor at the shock. 

3.2 The case of Bohm diffusion 

In this case we consider explicitly the most general situation in which the diffusion coeffi- 
cients and rate of losses are different upstream and downstream. The diffusion coefficient 
has the functional dependence Di^2{p) — -f^i,2P, where Ki^2 — (l/3)c^/(g-Bi_2) with c the 
speed of hght, q the charge of electrons and B the magnetic field in the shock region. 
This is also the same field strength used to determine the rate of losses: p — —Ai^2P^ — 
-(4/3)(TTf/B/(mec)V = -(4/9)^^f 2P^- We also have that Bi/ B2 = k<1. 



The function r reads: 



r,.fe/) = ^ln(l). (22) 

The general expression for the maximum momentum is then 
j(T-l\ ( 1 + Kr' \ ( 1 ~m(m.cf 

p_ = ^ (^ j j ( j ^jT^- (23) 

This value of the maximum momentum, though a reasonable estimate, is not to be 
interpreted as the exact momentum where a cutoff in the spectrum of accelerated particles 
should appear. From the physical point of view this can be expected since the residence 
times upstream and downstream, used to derive Pmaxi are defined as D/u^ only when losses 
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are unimportant, therefore by definition not at p ~ Pmax- IZirakashvih fc AharonianI (120071 ) 
found that for the case of Bohm diffusion the cutoff is located at 



namely 



Pr, 



3r 



-1/2 



3r 



For r = 4 and k, = 1 one has Pmax/po = 1-26. 



(24) 



(25) 



3.3 The case of Kolmogorov diffusion coefficient 

In this case D(p) = Kop^^^ and 



T{Z,Z 



2A 



(^2/3 _ ^/2/3) _ 



The maximum momentum is easily found to be: 

^3/4 

± Ui 17 — -L 

Pn 



SKoAr r + 1 



(26) 



(27) 



4 RESULTS 

In this section we discuss the results of the calculations. We apply the method to the three 
choices of the diffusion coefficient discussed above. We investigate the spatial and momentum 
dependence of the spectra of accelerated particles and the shape of the cutoff region for the 
three cases. 



4.1 Spectrum at the shock and shape of the cutoff 

In the left panel of Fig. [T] we show the spectrum of accelerated particles at the shock location 
for r = 4 and for the three choices of diffusion coefficient. The three diffusion coefficients 
are normalized in such a way to lead to the same maximum momentum Pmax- The curves 
refer to p'^Nq{p) and constant diffusion coefficient (dash-dotted line), Kolmogorov diffusion 
(dashed line) and Bohm diffusion (solid line). The spectrum of accelerated particles has a 
shallower cutoff in the case with constant diffusion coefficient. The cutoff becomes sharper for 
diffusion coefficients with more pronounced dependence on particle momentum. For instance 
the shape of the cutoff for Bohm diffusion is oc exp [—{p/po)'^], where Pq is defined in eq. [2H 
For the case of Bohm diffusion our resulting spectrum agrees well with the fitting formula 
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Figure 1. On the left (right) panel, spectrum of accelerated particles for r = 4 (r = 7) and for the three choices of the diffusion 
coefficient leading to the same value of pmax- The curves refer to constant diffusion coefficient (dash-dotted line), Kolmogorov 
diffusion (dashed line) and Bohm diffusion (solid line). The crosses in the left panel show the fit proposed by Zirakashvili & 
Aharonian (2007). 



proposed by IZirakashvili fc AharonianI (j2007l ). based on an interpolation of an analytical 



approximation valid for p > Pmax and numerical calculations in the transition region (shown 
as crosses in Fig. [1]). It is of some importance that the position of the cutoff is not exactly 
at Pmax but rather at po (see discussion in §3.21) . The discrepan cy between the calculated 



Zirakashvili fc Aharonian 



spectr a as presented here and as modeled by the fitting formula of 
( 120071 ) is typically fraction of percent, though at momenta p ~ Pmax it reaches ~ 10%, as 
visible in Fig. [H 

In the test particle limit the maximum compression factor which may be realized at the 
shock is r = 4, leading to Nq{p) oc p^"^. In cosmic ray modified shocks the effective total 
compression factor between upstream infinity and downstream can be much higher. For the 
sole purpose of illustrating the effect that may be expected in such cases, we consider here 
the case r = 7. The spectra, again normalized by the spectrum that would be obtained in the 
absence of synchrotron losses {No(p) oc in this case), are illustrated in the right panel 
of Fig. [T]for the three diffusion coefficients introduced above. At low momenta, p <^ Pmax, 
where synchrotron losses are negligible, the standard test particle spectrum is reproduced. 
The most noticeable difference with the case r = 4 is the well visible bump around p = Pmax 
for all choices of diffusion coefficient. The bump is broader for constant diffusion coefficient 
and becomes gradually more prominent and narrower for stronger dependence of the diffusion 
coefficient on momentum. 
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Figure 2. Spatial distribution of accelerated electrons for constant diffusion coefficient and r = 4. The left (right) panel refers 
to the upstream (downstream) region. The curves refer to p = Q.lpmax (solid line), p = pmax (dashed line) and p = 2pmax 
(dash-dotted line). 



4.2 Spatial distribution of the accelerated particles 

The approach described above also allows one to determine the spatial distribution of the 
electrons in the downstream and upstream regions, through Eqs. [7]and[Hl This information is 
important for the determination of the spectrum of the radiation emitted by the accelerated 
electrons. 

In Fig.[2]we show the spatial distribution of the electrons upstream (left panel) and down- 
stream (right panel) for constant diffusion coefficient and for three values of the momentum 
of particles: p = O.lpmax (solid line), p = Pmax (dashed line) and p = 2pmax (dash-dotted 
line). On the y-axis we plot p'^N{x,p). On the x-axis there is the distance from the shock 
in units of Li^max = Dq/ui upstream and L2,max = U2Tioss,2{pmax) downstream. In all cases 
a compression factor r = 4 has been adopted. 

It should be noticed that the diffusion length is independent of momentum in this case, 
therefore the size of the upstream region is basically independent of p at low momenta 
(compare the dashed and solid lines). At momenta close to the maximum momentum, losses 
become important, but by definition of maximum momentum the acceleration time and 
the loss time are comparable at that point, therefore the size of the diffusion region in the 
upstream fluid is roughly the same also for p ~ Pmax- At p > Pmax the density of particles is 
suppressed exponentially, as shown in Fig. [2] (dash-dotted line) . 

In the downstream region the situation is more interesting: there are two relevant spatial 
scales, = U2Tioss ~ 1/p and xfJJ = \/2D^tu^ ~ l/p^^"^. For a diffusion coefficient 
constant with momentum, the two spatial scales are equal at p^, = uI/{2ADq). Comparing 
this with the maximum momentum, one obtains: 
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Figure 3. Spatial distribution of accelerated electrons for a Bohm diffusion coefficient and r = 4. The left (right) panel refers 
to the upstream (downstream) region. The curves refer to p = 2pmax (dash-dotted line), p = Pmax (dashed line), p = O.lpmax 
(solid line). 
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Figure 4. Spatial distribution of accelerated electrons for a Bohm diffusion coefficient and r = 7. The left (right) panel refers 
to the upstream (downstream) region. The curves refer to p = 2pmax (dash-dotted line), p = pmax (dashed line), p = O.lpmax 
(solid line). The thin dotted lines show the result of a qualitative argument for the development of bumps, as discussed in the 
text. 



^ = -4^- (28) 

Pmax 2r(r-l) 

For r = 4, = 5/8, which imphes that at n < (5/8)pmax the size of the region where 

Pmax 

accelerated particles are located in the downstream fluid is of order xf^^^i^p) = U2Tioss{p)- This 
scaling with momentum can be qualitatively seen in the right panel of Fig. [21 by comparing 
the cases with p = O.lpmax and p = pmax- 

At small distances from the shock, p'^N{x,p) is spatially constant at low momenta, as it 
should be since synchrotron losses become negligible. 

Our results on the spatial distribution for the case of Bohm diffusion are illustrated in 
Fig. [31 The distances upstream are normalized to Li ^ax = D{pmax)/ui, while downstream 
L2,max = U2Tioss{pmax) ■ The left (right) panel refers to the spatial distribution of accelerated 



3 r + 1 
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particles {p'^N{x,p)) in the upstream (downstream) region. The curves refer to different 
momenta of accelerated particles: p = O.lpmax (solid line), p = Pmax (dashed line) and 
p = 2praax (dash-dotted line). 

The main differences with the previous case are all due to the fact that now the diffusion 
coefficient depends on momentum. In the upstream region, the size of the diffusion region 
for p < Pmax scales linearly with the momentum, as clearly shown in the left panel of Fig. [31 
At j9 > Pmax the synchrotron losses are faster than acceleration and the spatial size of the 
particle distribution starts decreasing. 

The situation downstream is again more interesting. Introducing once again the two 
distances xf^^^ = U2Tioss{p) ~ 1/p and xf^lf = \/1DqbPTioss{p) (constant in momentum), 
one sees that the spatial size of the downstream region for particles of momentum p is given 
by x1^'!^g = U2Tioss{p) oc 1/p for all momenta such that p/pmax < A/3(r + l)/(2r(r — 1)). For 
r = 4, the condition becomes p/pmax < 0.8. The scaling with 1/p can again be identified in 
the right panel of Fig. [3l Moving to lower and lower energies the drop in the spatial density 
of accelerated particles at x > L2^max becomes increasingly sharper. 

In Fig. m we plot the results for compression factor r = 7. The main difference with 
the previous cases is the appearance of pile-ups in the spatial distribution of low energy 
accelerated particles in the downstream plasma. This feature is due to the energy loss of 
particles at p ~ Pmax around the bump which appears right before the cutoff (see Fig. [1]). 
The bumps become more pronounced at lower energies while being absent at p > Pmax ■ 

The bumps are not artifacts of the calculation, as can be shown rather easily at least 
qualitatively: if advection is the only relevant process, then x = U2t downstream, and one 
can write dp/dx = —{A/u2)p^ ■ Conservation of the number of particles between the shock 
(where the particle density is iVo(p) and the location x, where the particle density is iV(p, x), 
then easily leads to 



Nq{p) p- 



■7 



N{p, x)dp = No{p)dp ^ N{p, x) = oc (29) 

( X _ ^ ) M _ 

where p = p/{^ — Axp/u2) is the momentum that a particle is produced with at the shock 
in order to reach the location x downstream with momentum p, as calculated from the rate 
of momentum losses. The last step in Eq. [29] holds for the case in which No{p) oc p~'^. In 
this case one can easily show that a spike develops at values of x such that the denominator 
vanishes, provided 7 < 2 (namely r > 4). In reality these spikes are smoothened by the 
presence of a cutoff in Nq{p) and even more important by the fact that at sufficiently high 
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Figure 5. Spectra of accelerated electrons upstream (left panel) and downstream (right panel) for a Bohm diffusion coefficient 
and r = 4. The curves are for x = —0.1L\^rnax (solid line), x = — Li.max (dashed line) and x = —2L\^rnax upstream and 
X = 0.1L2,maa; (solid line), X = L2^max (dashed line) and x = 2L2^max downstream. 



momenta diffusion prevails upon advection. In a qualitative way one could describe this effect 
by writing x = U2t + ^/2D{p)t and repeat the argument above. The result is not completely 
analytical but it requires the numerical integration of a simple differential equation which 
results in the dotted thin lines shown in the right panel of Fig. |H These qualitative results 
show how the bumps in fact are expected for r > 4 and for low enough momenta, while they 
disappear for p > p^ax- 



4.3 Spectra at different locations upstream and downstream 

In Fig. [5] we plot the spectrum of accelerated particles at different distances upstream (left 
panel) and downstream (right panel) for Bohm diffusion. In the left panel the spectrum is 
plotted at X = —O.lLi^max (solid line), x = —Li^max (dashed line) and x = —2Li^rnax (dash- 
dotted line), where Li^^ax = D{pmax) /ui. In the upstream plasma we can see the expected 
trend of lower energy particles to be confined closer to the shock surface. At increasingly 
larger distances from the shock only high energy particles are present. On the other hand at 
P > Pmax energy losses become important and the absolute number of particles drops down. 

In the right panel the spectrum is plotted at x = 0.1L2,max (solid line), x = L2,max (dashed 
line) and x = 2L2^rnax (dash-dotted line) with L2^rnax = U2Tioss{Pmax) ■ In the downstream 
plasma the main effect to be noticed is that the spectrum is truncated at lower energies 
while moving further downstream, as a result of synchrotron energy losses that cut out high 
energy particles from the spectrum. 
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Figure 6. Spectrum of accelerated particles for r = 4 and Bohm difFfusion, spatially integrated in the upstream fluid (dashed 
lines) and in the downstream fluid (solid lines). The left panel refers to the case in which no compression in the magnetic field 
is assumed, while the right panel shows the results for the case in which the magnetic field downstream is compressed by a 
factor vTT with respect to upstream. In both cases the momenta have been normalized to the same value of Pmax as calculated 
in the case without compression, in order to emphasize the fact that in the latter case the spectrum cuts off at a slightly lower 
momentum. 



4.4 Spatially integrated spectra of accelerated particles 



Although X-ray imaging has allowed us to observe the morphology of the radiation emitted 
by accelerated particles, it remains of interest to measure the spatially integrated emission 
from the shock region, which often is the only quantity we can actually measure. In this 
section we discuss our results for the spatially integrated spectra of accelerated particles, 
which have a behavior similar to that of the integrated emission. The spatially integrated 
spectra can be obtained by integrating Eqs. [7]and[H]in x in the appropriate regions upstream 
and downstream, namely: 



Slip) 



dxNi {x, p) 



S2{P) 



dxN2{x,p). 



After calculating analytically the integral over x we obtain: 



Slip) 



2^7 



dp'NQ{p')\l + Erf 



dp'D^{p')No{p') X ;^exp 



AAir 



Ui[Z — z) 
2Airi/2 



<Pi{.p') 



Proceeding in the same way, the integral over the downstream region reads: 



S2{P) 



1 



dp'D2{p')No{p') X ;^exp 



uljz'-z)' 



z'^U2 

2l7 



dp'No{p') 



Erf 



U2{Z - Z) 

2A2ri/2 



Mp')- 



(30) 



(31) 



(32) 



Here Erf{y) is the error function. 
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These expressions apply to any choice of the momentum dependence of the diffusion 
coefficient. The quahtative features of the integrated spectra upstream and downstream for 
V < Vmax can easily be understood. Let us for instance consider the case of Bohm diffusion: 
the upstream size of the region where particles with momentum p diffuse is ~ D{p)/ui oc p. 
Therefore Si{p) oc Nq{p)p oc p~^ for r = 4. The spatial size of the region where particles 
reside in the downstream region is ~ U2Tioss,2 oc 1/p, therefore S2{p) oc Nq{p)u2Tioss,2 oc p~^ 
for r = 4. Again, this trend is easily seen in the results of our exact calculations, shown 
in the left panel of Fig. [6l The solid line shows the spectrum as spatially integrated in 
the downstream region, while the upstream integrated spectrum is shown as a dashed line. 
The sum of the two is the dash-dotted line. Since Si{p) is much flatter, at most energies 
the largest contribution to the spectrum of radiating electrons in an astrophysical source 
comes from the downstream region. A small contribution appears only at p ~ Pmax- This 
effect is even more visible if one assumes that the magnetic field downstream is obtained 
from the compression of the magnetic field in front of the shock surface, by an amount k~^. 
In the right panel of Fig. [6] we show the case in which = vTT, which is the case of 
compression of an isotropic turbulent field upstream for r = 4. One can see that in this case 
the downstream component produces a more pronounced bump in the total spectrum, as 
due to the fact that particles immediately downstream at p ~ Pmax lose energy more rapidly 
than upstream. While being hardly visible in the spectrum of synchrotron emission because 
this is dominated by the downstream region, this feature might be visible in the spectrum 
of inverse Compton radiation in those cases in which gamma ray emission can clearly be 
attributed to this process rather than to a hadronic origin. 

The simple trends illustrated above can be easily extended to any dependence of the 
diffusion coefficient on momentum. For instance for a constant diffusion coefficient, one can 
expect that Si{p) oc p^"^ and S2{p) oc p~^. The slope of S2{p) is independent of the diffusion 
coefficient as it is determined by energy losses at all energies. Of course this simple prediction 
is the consequence of the assumption of stationarity of the problem. In real astrophysical 
sources, stationarity will be reached only at sufficiently high energies, while at lower energies 
the size of the downstream region is determined by the age of the source and is therefore 



independent of momentum. This is t 
spectra of radiation (see for instance 



le reason for the appearance of spectral breaks in the 



KardashevI (119621 )). 



Another note of caution should be issued about the assumption of planar shock adopted 
throughout the calculations: in astrophysical sources the effects of sphericity might become 
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important, especially on the volume integrated quantities. On the other hand one may argue 
that since particles only move a distance ~ U2Tio.ss,2{E) ocl/E downstream before losing an 
appreciable fraction of their energy, the effects induced by spherical symmetry are expected 
to become more important at low energies. It is likely however that the assumption of 
stationarity breaks down before the spherical geometry leads to important effects, but this 
should be checked case by case. 

5 CONCLUSIONS 

We proposed a semi-analytical calculation of the process of diffusive acceleration of electrons 
at a non-relativistic shock in the presence of synchrotron (or inverse Compton scattering) 
losses. The calculation returns the spectrum of electrons at the shock and at any location 
upstream and downstream for an arbitrary choice of the momentum dependence of the 
diffusion coefficient. The results of the calculations for the spectrum at the shock have 
been illustrated for three cases: 1) diffusion constant in momentum {D{p) = Dq), 2) Bohm 
diffusion {D(p) oc p), and 3) Kolmogorov diffusion {D{p) oc p^^^). For the space dependence 
of the spectrum and the integrated electron spectra we restricted for simplicity to cases 1) 
and 2). 

While confirming the results of previous calculations for the simple case of a constant 
diffusion coefficient, our results can be easily obtained for any choice of the diffusion coeffi- 
cient. Most results illustrated here are referred to the case of Bohm diffusion and to different 
compression factors at the shock. While the slope at low momenta is the standard test par- 
ticle slope, depending only upon the compression factor, the shape of the cutoff, which is 
determined by the onset of synchrotron losses, depends on the adopted diffusion coefficient. 
It is exponential (oc exp [—p/po]) for D{p) = Dq and approaches a oc exp [—{p/po)^] for 
Bohm diffusion, where po is related to Pmax through Eq. [25l 

The importance of these calculations for the description of the phenomenology of su- 
pernova remnants, and possibly other classes of sources, is evident: for supernova remnants 
the synchrotron X-ray emission is now resolved both spectrally and spatially and its care- 
ful description could allow to access precious information on the acceleration process. For 
instance, if the magnetic field is indeed amplified by accelerated particles (mainly protons) 
at the shock, the maximum energy of the electron component is determined by synchrotron 
losses, while if no amplification occurs, the maximum energy could be determined by the 
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finite age (or spatial size) of tlie accelerator, as it is for protons. The two cases result in 
different shapes of the cutoff in the electron spectrum (exponential in the first case and 
super-exponential in the second case) and possibly in different fits to the synchrotron X-ray 
flux and morphology. 

One note of caution should be issued about the assumption of stationarity that is under- 
ling all calculations of this type. Strictly speaking, stationarity can be reached only if the 
time for energy losses is shorter than the age of the shock at all energies and the maximum 
energy is roughly obtained by equating the acceleration time with the loss time. In a real 
astrophysical source it is usually the case that at sufficiently low momenta the loss time is 
shorter than the age of the source, so that no stationarity can be actually reached at those 
momenta. In other words a spectral break can be expected in the volume integrated spectra 
of accelerated particles at some momenturn Vh an d the stationary solution found here should 
be applied only for p > pb (e.g. (lKardashevlll962l )). This reflects in the spectra of the emitted 
radiation. 

The formalism presented in this paper will be extended to the case of cosmic ray modified 
shocks in an upcoming paper: this application is crucial in that the presence of a precursor 
may appreciably flatten the electron spectra at high energy and lead to the production of 
more pronounced spectral bumps. 



ACKNOWLEDGMENTS 

The author is grateful to L. O'C. Drury for reading a preliminary version of the manuscript 
and E. Amato for a useful conversation on the spatial distribution of accelerated particles. 
This work was partially supported by MIUR (under grant PRIN-2006) and by ASI through 
contract ASI-INAF 1/088/06/0. This research was also supported in part by the National 
Science Foundation under Grant No. PHY05-51164, in the context of the Program Particle 
Acceleration in Astrophysical Plasmas, July 26-October 3, 2009 help at the KITP in Santa 
Barbara. 

APPENDIX A: GREEN FUNCTION OF THE ADJOINT EQUATION 

Here we briefly illustrate the determination of the Green function of the adjoint equation, 
defined by: 
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dG (9^^ dG 

where Q = Q{x,p;x',p') and A and -u are the appropriate quantities for the upstream and 
downstream regions. We introduce the following variables: 

and 

1 r, oiv) 



In terms of these new variables, the equation above becomes: 

de dr Ap'^ ' ^ ' 

Fourier transforming this equation with respect to ^ one easily finds the solution in the form: 
1 



^(0 



\ 1/2 


" e' 


) exp 


4t_ 



(A5) 
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